Superintegrable quantum u (3)— systems and 



higher rank factorizations 

J. A. Calzada*, J. Negro*, M.A. del Olmo* 

t Departamento de Matemdtica Aplicada 
* Departamento de Fisica Teorica, Atomica y Optica 
Universidad de Valladolid, E-47005 Valladolid, Spain 

E. mail: juacal@eis.uva.es, jnegro@fta.uva.es, olmo@fta.uva.es 

February 7, 2008 
Abstract 

A class of two-dimensional superintegrable systems on a constant curvature surface is 
considered as the natural generalization of some well known one-dimensional factorized 
systems. By using standard methods to find the shape-invariant intertwining operators 
we arrive at a so(6) dynamical algebra and its Hamiltonian hierarchies. We pay attention 
to those associated to certain unitary irreducible representations that can be displayed by 
means of three-dimensional polyhedral lattices. We also discuss the role of superpotentials 
in this new context. 

1 Introduction 

This work deals with a class of superintegrable Hamiltonians systems, in the framework of 
the Schrodinger equation of quantum mechanics, and its connections with the factorization 
method. We will restrict ourselves to a particular case where the underlying symmetry is the 
Lie algebra u(3), but its main features can be directly implemented to higher dimensional 
systems. 

The main objective of this study is to show a natural extension to higher dimensional 
spaces of the intertwining (or Darboux) transformations from a well known class of one-dimen- 
sional factorized systems. In fact, we want to set the higher rank «(n)-systems corresponding 
to those having as dynamical algebra the Lie algebra of rank one u(2). We will show in detail 
that the application of procedures familiar in one dimension to a concrete two-dimensional 
system will lead us to a wide set of operators closing a dynamical Lie algebra. We also consider 
discrete symmetry operators quite important to perform equivalences. All these operators 
connect eigenstates that can be drawn as points in a three-dimensional lattice giving rise to 
polyhedrons representing degenerate series of u(3) irreducible representations. Each of these 
series corresponds to the same energy and can be embedded in just one representation of the 
Lie algebra so(6). 

The notion of superpotential will also be re-examined inside the higher rank formalism. 
Thus, the usual procedure to look for solutions with separable variables can be better appre- 
ciated under this point of view. 



Thus, we try to implement the program of generalization of the factorizable one-dimensional 
systems involving Lie algebras of rank one as dynamical algebras, as can be seen, for instance, 
in the classical paper by Infeld and Hull pQ. We also hope that this work will be useful when 
dealing with other integrable systems, but not necessarily maximally integrable, for instance 
not enjoing for such a wealth of factorizations, or even not having a system of separable 
variables, but still allowing for algebraic methods jH El EJ 03 El d • 

The organization of the paper is as follows. In section [21 we will introduce a two- 
dimensional superintegrable system and find some separable solutions by standard proce- 
dures. Although these polynomial solutions are known and can be found in other references, 
this will serve us to recall some aspects of the usual factorization technique and to precise 
the operators of the Lie algebra u(2) related with the Lie algebra so(4), and their support 
spaces. Next, in section|31 we will look for other sets of intertwining operators, corresponding 
to the Lie algebras u(3) and so(6), taking also into account discrete symmetries. We char- 
acterize the eigenfunctions belonging to irreducible representations that will be depicted as 
the points on octahedrons and the interpretation of some of its planar sections. The analog 
of the superpotentials and their relation with certain types of solutions will be considered in 
section^ Some conclusions and perspective for future work will close the paper. 

2 A superintegrable w(3)— Hamiltonian system 

We will fix our attention on a superintegrable Hamiltonian system defined inside a three 
dimensional Euclidean ambient space (SH^IEH^D- In fact, our system lives on the 2-sphere 
S 

S = (so) 2 + (si) 2 + (s 2 ) 2 = 1, (s , si, s 2 ) € M 3 
In the frame of the Schrodinger equation, this Hamiltonian takes the form 

V 2 J (s ) 2 (si) 2 (s 2 ) 2 

where (lo,h,fa) G M 3 , and Ji = —eijkSjdk (note that the Jj's operators generate the rotation 
Lie algebra so(3)). We can parametrize S by means of spherical coordinates (fa, fa) around 
the s 2 axis given by 

so = cos fa cos 0i, si = cos fa sin fa, s 2 = sin fa (2-2) 

Then, the eigenvalue problem 

H<$> = E$> 

after substituting the coordinates (|2.2|) . takes the form of a separable differential equation 

<$> = E<$> (2.3) 
The solutions separated in the variables fa and fa, i.e. 

*(0i,02) = f{fa)g{fa) 

2 



I 2 — 1/4 

-9| a + tan fadfa + — 



+ 



1 



sin fa cos 2 



cos 2 fa 



-1/4 



sin 2 fa 



after replacing in (|2.3|) originate the equations 
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where a is a separating constant. Next we will solve each of these two equations through 
standard factorizations giving rise to polynomials. The key point is that the results obtained 
for the first equation will match in a certain way with those of the second one originating 
degenerate levels. 



2.1 The fa— factorization 

The one-dimensional Hamiltonian (J2.4|) in the variable <p\ is a well known example in the 
theory of factorizations pQ. So, in the following we will restrict ourselves to give a list of the 
relevant results. We will see later, in section |SJ how to make use of these considerations in a 
broader context. 

The second order operator at the l.h.s. of eq. (|2.4|) can be cast as a product of first order 
operators 

being 

Aq = ±8^ - (Z + 1/2) tan fa + (h + 1/2) cot fa A = (Z + h + I) 2 
These elements are part of a family of operators {i^, A~ , A m , Hf\}, m € Z, where 

A± = ±8^ - (Z Q + m + 1/2) tan fa + (h + m + 1/2) cot 0i (2.6) 
A m = (Z -Hi + 2m + l) 2 
01 2 (Z + m) 2 -l/4 (Z 1+ m) 2 -l/4 

They originate the one-dimensional Hamiltonian hierarchy (J2.7JI . starting from #(oV r ^ ne 
Hamiltonians satisfy the fundamental relation 

^(.m) = A m + A m = A m _ l A^ n _ l + A m _i (2-8) 
so that are shape invariant intertwining operators, i.e. 

An-H(m) = ' ^ro#(m+l) = H {m) A m ( 2 - 9 ) 

Hence, from a formal point of view, the operators acting on a Hamiltonian eigenfunction 
will give another eigenfunction of a consecutive Hamiltonian in the hierarchy with the same 
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eigenvalue. If we design the eigenfunction spaces of Ht\ (as differential operators) by Tif^ , 
then we have 

A— . n-(<t>i __, n-/4>i a+ . -ij^i -yi<t>\ 



In principle, the discrete spectrum and the physical eigenstates of H?q\ could be obtained from 
the fundamental states /P» and their eigenvalues of all the Hamiltonians in the hierarchy 

{Ht\y. These fundamental states are determined by /?> = 0, giving the solutions (up 
to a normalization constant) 



/ ( ° m) (0i) = cos' 0+m+1 / 2 0i sm'^/a , 



with eigenvalues A m = (Zo + h + 2 m + l) 2 . Often the intertwining operators Q2.6|) are written 
in the form 

At = ±6^ + u m (<h), "m(&) = ^ 1 (2.10) 

where w m (0i) is called superpotential function. 

In order to go from the ground eigenstate, fP n of Hf 1 *., up to the excited eigenfunction, 
of , with the same eigenvalue, we apply consecutive operators A + 

/(o) = A$ A£ ■■■ A^-x /( m ) ( 2 -H) 

obtaining explicitly 

fffa = N sm h+1 / 2 fa cos' 0+1/2 ^i P,g 1 ' ,o J(cos(2 0i)) (2.12) 

where Pi a ' b \x) are Jacobi polynomials and iV a normalization constant. Therefore, the 
spectrum of the first separating Hamiltonian (|2.4|) is given by 

a = \ m = (Z + in + 2 m + l) 2 , meZ + (2.13) 

The following two subsections are devoted to characterise the Lie algebras of shape in- 
variant intertwining operators for the one-dimensional Hamiltonian hierarchies. They will 
constitute a useful pattern for the two-dimensional Hamiltonians of section |3J 

2.2 The dynamical algebra u(2) 

Starting from the operators A^ let us define free-index operators A^ acting inside the total 
space ® m 7^ m , in the following way [T2l IT3*] : 

A + /(m+l) := 2^m/(m+l) f(m) 

A~f(m) ■= 2^m/(m) a /(m+l) ( 2 -14) 
^/(m) : = -^o+ii+2m)/ (m) a / (m ) 
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where /( m ) (or f( m )) denotes an eigenfunction of Hr m y This action can be extended to linear 
combinations of eigenfunctions by linearity. With this convention we can rewrite (|2,8|) and 
(|2.14j) simply as the commutators 

[A, A ± ] = ±A ± , [A-,A + ] = -2A (2.15) 

assuming that the action is on any (linear combination of) /( m )'s. The commutators ()2.15j) 
close the Lie algebra su(2), whose Casimir element is given by C = A + A~ + A(A — 1). The 
eigenvalues of C, labeling the irreducible unitary representations (IUR), are j(j + 1), where 
2j G Z + . The dimension of the support spaces of these IUR's is, obviously, 2j + 1. We 
make use of the standard notation | j, s) for an ^-eigenvector with eigenvalue s, inside the 
'j-representation'. 

Now, we can identify the eigenstates of the Hamiltonians Ht\ in terms of representation 
vectors |j, s). First, let us consider the ground states /? > characterized by 

A-f {m) = 0, Af° m) = -[(l + + 2m)/2]/ ( ° m) (2.16) 

following the notation (|2.14|) . These relations suggest the identification (up to a normalization 
constant) 

/(m) = \3m, -jm), jm = (k + h + 2m)/ 2 

To see that indeed this is the case we need to define the whole representation space as well 
as an inner product. Thus, consider the space £ 2 [0, 7r/2] of square integrable functions in the 
interval [0, 7r/2]. Then, the wavefunctions obtained from the ground state /L by the con- 
secutive action of the operator A + will span the representation space of a j m -representation, 
with j m = (Iq + l\ + 2m)/2, provided that both Iq + m and l\ + m belong to Z + . The 
wavefunctions of the space so generated vanish at the end points (i.e. ^(0) = ^(tt/2) = 0), 
and the hermiticity relations (A~)^ = A + , A^ = A are implemented in all the space. Hence, 
under these conditions, (A + ) k f® m ^ can be identified, up to normalization, with the vector 
state |jm> jm 

+ k). 

As a consequence, the excited states obtained in this way for any Hamiltonian in a fac- 
torization hierarchy where Iq and l\ are positive integers, correspond to IUR-vector states. 
For instance, the eigenstate of the fc-th excited level of is 

f( ) = \jk + k,-j k + k), j k = {l + h + 2k)/2, k = 0,1,2... 

and #/q\ (as well as any Hf\) can be expressed in terms of the sTi(2)-Casimir C acting on 
such representations 

Hf Q ) = 4(C + 1/4) 

Therefore, the eigenvalue equation for any of the excited states can be written as follows 

ffgj/fo = A(C + l/4)\j + k,-j + k) 

= 4(j + k + l/2) 2 |j + k, -j + k) = (l + h + 2k + l) 2 /f ) 

with ft = 0,1,2,.... 
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It will be convenient to consider a new diagonal operator D, to be added to the generators 
of su(2) ET1 . define by 

Df( m ) '■= Go - h)f(m) 

It is immediate to see that D commutes with any other operator of su(2) giving rise to the Lie 
algebra u{2). In this way any eigenstate in the Hamiltonian hierarchy can be characterised 
completely by an eigenfunction of an n(2)-IUR. Without D we would have an ambiguity 
due to the fact that different fundamental states with values of Iq and l\ giving the same 
jo = (Iq + h)/2 would lead to the same ^-representation of su(2). 

It is worthy noting that when 1$ or l\ are not in Z + the eigenfunctions and spectrum of 
the Hamiltonian hierarchies are still given by (|2.12j) and (|2,13|) . but, these states belong to 
non-unitary representations of u(2). 

2.3 The dynamical algebra so(4) 

As we have just seen in the previous subsection the eigenstates sharing the same energy of 
the one-dimensional Hamiltonian hierarchies in the variable </>i are given in terms of IUR's of 
the dynamical algebra u(2). However, in this respect, there is a point not quite satisfactory: 
different u(2)-IUR's may correspond to states with the same energy. We would prefer a larger 
dynamical algebra with a simpler correspondence, i.e., such that only one of its IUR's gives 
all the eigenstates with the same energy in the hierarchy. 

In order to build up a dynamical algebra having these properties, let us introduce the 
two-dimensional parameter space (Zo, h)- Any operator with one subindex defined in subsec- 
tions 12.11 and 12.21 will change to a two-subindex notation in the following way: 



1. The one-dimensional Hamiltonian (|2.4j) will be denoted by Hn Q 



(io,h) 



Its eigenfunctions will be designed by /a 0) M- 

The factor operators in (|2.6j) will be rewritten as A^ o j » 

Af, M = ±9^ - (Z + 1/2) tan 0i + (h + 1/2) cot^, An,, = --(l +h] 



Go.ii) ^ ~ 10 ^ 1 ' Lmi ^i i" W ™<Pi, ^(/ ,ii) - 

Now, in this way, relations (|2.9|) can be expressed as 

a— u4>i _ tj4>i a— a+ u<t>i _ jj4 i i a+ 1 7\ 

^(io.ii) (io.ii) _ (io+l.ii+^^io.ii) ' ~ n {loM {lo,h) 1 '> 

With this convention we can also define the free-subindex operators A^ 1 , A, D as in 1)2. 14|) . 

On the other hand, notice that each two-parameter Hamiltonian Hfo j\ is invariant under 
the reflections 

: (hih) —* (—lo,h), h ■ {lo,h) — » (lo,—h) 



G 



This property gives rise to a second factorisation (see also |2JE3 El) via conjugation of the 
operators of the first factorisation by the reflection opera tors 

IqA^Iq = i± I AI = A, IqDIq = D 

hA^i = i T , hAh = -A, hDh = -D 

Explicitly 

Af kM = ±8^ + (lo-l/2) tanfc + (h + 1/2) cot^, = ~(-l +h) (2.18) 

The above operators {^4, A ± } generate a Lie algebra isomorphic to su(2) denoted by su(2). 
Since su(2) and su(2) commute and, essentially, D and D coincide with A and A, respectively, 
the complete dynamical algebra has the structure of a direct sum su(2) © su(2) rs so(4). 

If we allow to act with the so(4) generators on an Hamiltonian Ha ^\ we will get a two- 
dimensional parameter lattice of Hamiltonians which constitute a so(4)-hierarchy fixed by 
the initial values (Iq, {Hi _ n+m ^ 1+n+m \ , m,n £ Z. Each energy level of this Hamiltonian 
hierarchy is degenerated and the eigenstates belong to so(4)-representations. 

Let us concentrate on the hierarchies associated to IUR's of so(4). Now, these so(4)-IUR's 
are fixed by the fundamental (or lowest weight) states satisfying 

A Ilo M fioA) = A Ii oM fioM = ( 2 - 19 ) 
These are realized, up to a constant, by the wavefunctions 

/ ( ° 0in) = cos 1 / 2 0! sin" +1 / 2 fa, n G Z + (2.20) 

where we have taken Iq = and l\ = n. We see also that the state (|2.2U|) is stable under 
Jo (i.e. Iofn = fn jjO) an d comes into the other fundamental state (annihilated by A + 

and A + : the highest weight) of the same representation. Hence, these representations will 
be invariant under Iq and I\. Therefore, the so(4)-IUR's obtained from Q2.2UJI are symmetric 
tensor products that can be denoted by 

j = h/2 = n/2, nGZ^° 

where l j } stands for a ^-representation of su(2). In this way the degenerancy of the re-th 
energy level is (n + 1) x (n + 1), which is composed of n + 1 IUR's of u{2) each of them of 
dimension n + 1. 

The Hamiltonians in this hierarchy can be expressed in terms of any of the s«(2)(or 
su{2)) Casimir operators H^ lo>l ^ = 4(C + 1/4) = 4(C + 1/4). With the help of all the discrete 
reflections we get directly its expression also in terms of the so(4)-Casimir 

H {khh) =(C + 1/4) + I (C + l/4)/ + h(C + 1/4)/! + I h(C + l/4)/ /i 

= {A + ,A-} + 2A 2 + {A + ,A-} + 2A 2 + l 

= {A + ,A-} + {A+,A~} + L 2 + L 1 2 + l 
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where the diagonal operators L and L\ are defined by 



Certainly, some so(4)-hierarchies (those corresponding to the IUR's previously described) 
may have Hamiltonians whose explicit expressions coincide 



and the same happens with their corresponding eigenstates. But we can not get rid of this 
multiplicity unless we enlarge the ambient space. 

Another natural question is whether there are other intertwining shape- invariant operators 
inside the so(4)-hierarchy. We can build, for instance, other pairs of operators through the 
composition of those already known 



This kind of shape-invariant operators change two units either the parameter Iq or l\ (but 
not both at the same time). When we restrict to Iq = or l\ = there are also first order 
intertwining operators changing one unit the nonvanishing parameter. This feature is not so 
special; it is also shared by the 'radial oscilator' hierarchies (which are closely related to 
the ones presented here). 

For other Hamiltonian so(4)-hierarchies the physical eigenstates are described by non- 
unitary representations that are not invariant under both reflections. In this respect, their 
description becomes more involved, so that one must be very careful in these cases. 

2.4 The 2 — factorization 

Now, let us return to the separation process started in subsection 12.11 The second equation 
(|2.5|) obtained from the initial separation of variables can be dealt with along the same lines, 
substituting the eigenvalues obtained from the previous factorization, a = \ m = (Iq + l\ + 
2m) 2 . The most relevant fact, here, is that the new factorization leads to a degeneration of 
the energy levels which suggest that the underlying dynamical symmetry could be larger, as 
it will be confirmed in the next section. Thus, substituting in (|2.5j) . we have 



H (hM - H (-l +l,h) ~ H (ku-h+l) 



X ± = A ± A ± 



y 





{d H - (k + h + 2(m + 1)) tan(0 2 ) + (Z 2 + 1/2) cot(0 2 )} 
x {-3<fe - (h + h + 2m + 1) tan(0 2 ) + (h + 1/2) cot(<£ 2 )} 
+ {h + h + h + 2m + 3/2) (Z 2 + (l + h + 2m + 5/2) 



(2.21) 



M+M ~ + Mo - 




in the variable c/> 2 



1 + V-n-l 
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where 

M± = ±d H - (l + h + 2{m + 1) + n) tan(0 2 ) + (h + n + 1/2) cot(0 2 ) 
= (h + l + h + 2n + 2m + 3/2)(Z 2 + ii + io + 2n + 2m + 5/2) 

Now, the values for the energy (following closely the same arguments of section |2~T|) are given 
by 

E = fi n = (l x +l + l 2 + 2n + 2m + 3/2)(Z 2 + h + l + 2n + 2m + 5/2) (2.22) 
The fundamental states g® n ^ for this factorization are 

5 o n) (0 2 ) = Nco6^ +l0 ^ +2m+1 <h sin' 2+n+1 /2 2 

and the eigenfunctions of the initial Hamiltonian (|2.21f) can be written in the form 

5 ™ o) (0 2 ) = cos h+l0+2m+1 4> 2 sin i2+1 / 2 4> 2 Pi l2+1 ^ h+k > +2m+1 \cos2 ( p 2 ). (2.23) 

The commutation relation for the relevant free-index operators M 1 * 1 , defined in a similar way 
as A ± in Q2.14JI . is again that of su{2), 

[M~,M + ] = — 4(Zi + l + l 2 + 2m + 2n + 1) = -2M 

The eigenfunctions (|2.23f) are square- integr able, but the representations are unitary provided 
that, besides the previous conditions on Iq and h, the parameter l 2 be also a positive integer 
number. 

In summary, if we finally join the results of both factorizations, the square-integrable 
eigenfunctions of the Hamiltonian (|2.3|) in the separable variables (4>i,4> 2 ) are given by the 
products 

<&m,n(0i,02) = /(^)(0i)5(o) (<h), m,neZ + (2.24) 

where the components have the polynomial expressions (|2,12|) and (|2.23|) . The corresponding 
eigenvalues given in (|2.22|) are degenerated for those values of m and n whose sum m + n 
keeps constant (see also, for instance, Ref. jllj). 

3 Dynamical symmetries 

The spectrum obtained by the methods of section |^1 suggest the existence of a bigger dy- 
namical algebra of the Hamiltonian hierarchy. This is the point that we want to address 
here developing exhaustively the concept of intertwining (shape invariant) operators for this 
kind of Hamiltonians. Such operators will supply us with a more consistent picture of the 
spectrum and eigenfunctions. Thus, based on the considerations of subsections 12.21 and 12.31 
we will introduce three sets of intertwining operators closing the Lie algebra u(3). Then, 
in the following subsection, we will enlarge this algebra to so(6) by means of the relevant 
reflections. 
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3.1 The Hamiltonian u(3)— hierarchies 
3.1.1 The set {A+, A , A} 

As we will use some properties of section [2] in a different direction, it is convenient to intro- 
duce another notation more appropriate to rewrite some previous results. The Hamiltonian 
(|2.1|) characterized by the parameters I = (Zo, ^1,^2) will be referred to as Hr^^i^, and the 
operators defined by (|2,6j) will be taken henceforth with a three-fold subindex 

Af lo>hih) = ±6^ - (Z + 1/2) tan fa + (h + 1/2) cot fa (3.1) 

Since the differential operators (|2.6|) and <|3. 1|) depend only on the variable fa, they do not 
affect the part in the total Hamiltonian (|2.1|) depending on the second separable variable fa. 
So that, in the same way as Q2.17JI we have the intertwining relations 

A (lo,h,h) H ( l o,h,h) = H (lo+hh+l,l2) A (l ,h,l2) 

This means that now A7 1q ^ ^ is acting on eigenstates of H^^^ leading to eigenstates 
of fl'(j[ +i,j 1 +i ) j 2 ), while A^ o ^ ^ does it in the opposite way (later we will comment on the 
square- integrability conditions through unitary representations). 

If we include the normalizing constant just as in (|2.14| . and define global operators acting 
on eigenfunctions of this class of Hamiltonians in the form 

^ + *(io+l,ii+i,J2) := 5^fo,Ji,j a )*(fo+:Mi+l,la) x *(fe,ii,fa) 

A ®(lo,h,h) := + 'l)*(Wi,fe) 

we are lead to the standard su(2) commutators (|2.15() . Here, we want to stress again that now 
these operators are acting on the total wavefunction of complete Hamiltonians like Hn Qi i lt i 2 \, 
not just on a factor function in only one variable. 

In order to introduce other sets of operators we will use the fact that the Hamiltonian l|2.1j) 
can be separated in other coordinate systems. Since the axes (so, si, S2) P^ a y a symmetric 
role in the Hamiltonian, we will take their cyclic rotations to get two other sets of coordinates 
and, hence, new sets of intertwining operators. 



3.1.2 The set {B+,B-,B} 

We will take the spherical coordinates choosing as third axis not S2, but s±, i.e. 

S2 = cos £2 cos£i, so = cos £2 sin£i, si = sin £2 (3-2) 

Then, the initial Hamiltonian is also separated in the coordinates (^1,^2)- I n particular, we 
can build the operators ^ ^ in a similar way as A^ ^ ; %. From the coordinate systems 
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'2 

(3.3) 



(|2.2|) and (|3.2|) we easily arrive at the following expressions for the new set in terms of the 
initial coordinates (0i,<^>2) 

B^ o k ^ = ±(sin <f>i tan^>2^i + cos ^l^) ~~ (^2+1/2) cos ^icot^ + (Zo+1/2) sec^itan^; 

These operators intertwin the pair of Hamiltonians 

B {l ,h,h) H ( l o^M = H (lo+l,h,h+l) B {l ,h,l 2 ) 
B (l ,h,h) H (lo+l,h,h+l) = H (lo,h,h) B (l ,h,h) 
The 'global' operators, defined by 

B+ ®(l +l,h,h+l) := 5 B (t ,W 2 ) # (fo+l,M 2 +l) « *(Wi.fa) 

B ~*(Wi.*a) := h B (l ,li,h)®( l °' l iM K *Po+l.'i.'fl+l) 
B ®(lo,h,h) := ~Wo + ^)*(ZoA,Z 2 ) 
also close a new su(2). 

3.1.3 The set {C+,CT,C} 

Finally, taking the spherical coordinates around the so axis, 

s\ = cos 6*2 cos 0i, S2 = cos #2 sin#i, so = sin 02 

the Hamiltonian is also separated in the variables {61,62} and we get a new pair of operators, 
that written in terms of the initial <j>i and <fo variable, take the expression 

Cq q1i l2 \ = ±(cos </>i tan 02^0! — sin</>i<90 2 ) + (Zi— 1/2) cosec0itan</>2 + (Z2+I/2) sin0icot</>2 

(3.4) 

These operators act as intertwiners of the Hamiltonians in the following way 

C (lo,h,h) H 0o>h-l,h+l) = H (lo,h,h) C (l ,h,l2) 
The 'global' operators are defined by 

closing the third algebra su(2). Notice that C = S — A. 

In fact, as we saw in sectional each separable system gives rise to two sets of intertwining 
operators (in that section distinguised by means of the tilde). However, here we have made a 
'good' choice of the above three sets that will close a Lie algebra (on this point see section l3~2|) . 
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3.1.4 The complete algebra u(3) 



Now, we can join all the transformations above defined, A^, A, B^,B, C^, C, and commute 
any two of them to check that indeed they close a Lie algebra su(3). The nonvanishing 
commutators are 

[A 3 ,A ± ] = ±A ± [A-,A+]=2A [A+,B-] = C- [A+,B} = -A+/2 
[A+,C+] = -B+ [A + ,C]=A + /2 [A-,B + ] = -C + [A-,B]=A~/2 
[A + ,C + ]=B- [A~,C] = -A-/2 [A,B + ]=B + /2 [A, C + ] = -C + /2 
[A,C-] = C~/2 [B,B ± } = ±B ± [B-,B + ]=2B [B+,C-] = -A+ 
[B + ,C} = -B+/2 [B-,C+} = C+/2 [B-,C] = B-/2 [B,C + ] = C + /2 
[B,C-] = C-/2 [C7,C ± ] = ±C ± [C + ,C-] = 2C [A~,C~} = B~ 

The Casimir operator is given by 

C = A + A- + B + B- + C + C- + -A{A - 3/2) + -B(B - 3/2) + \c{C - 3/2) (3.5) 

o o o 

In order to complete an algebra n(3) we can add a diagonal operator D commuting with all 
the above transformations. It is a central operator, i.e. 

D:=l -h-h, [D,-]=0 

We can also adopt the global operator convention H for the Hamiltonians in the hierarchy 
by defining its action on the eigenfunctions j 2) / 3 ) of ^(j 1) / 8l z 3 ) by 

In this way we can express the Hamiltonian H in terms of both operators C and D 

1 15 

H = AC--D 2 + - (3.6) 

In the case of one-dimensional systems, one (first order) intertwining set {^^j for the Hamil- 
tonian gives rise to its factorization. However, for Hamiltonians with more degrees of freedom 
(more components, or in more dimensions) the relationship of H with these operators, in 
general, turns out to be more complex. In our case the set {^4 ± , B^ , C ± } according to ex- 
pressions p.5(l and ()3.6|) is enough to express the Hamiltonian as a certain quadratic function 
H = h(A + A~, B + B~ , C + C~) generalizing the usual factorization. 

In summary, we have built an algebra u(3) of intertwining operators that, once fixed the 
initial Hamiltonian with parameter values (lo,h,h), gives rise to a two-parameter Hamilto- 
nian hierarchy 

{^(lo+m,h+m-n,l2+n)}j 171, n £ 1, 

where the points (Iq + m, l\ + m — n, I2 + n) lie on a certain plane D = d,Q. In this subsection 
we will consider this special hierarchy, together with its eigenstates, connected to the IUR's 
of u{3). The states of such representations are square integrable and, therefore, should take 
part of the physical eigenfunctions whose energy eingenvalues belong to the spectrum. 
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In order to build an IUR we start from a fundamental state <3? annihilated by A and C 
(two simple roots of su{3)) 

Aj$ t = Ci$ £ = (3.7) 
with £ = (lo,h,h)- Such states exist only when l\ = 0, taking the explicit form 

$l(<h,<h) = N cos lo+1/2 fa sin 1/2 0i cos' 0+1 cp 2 sin' 2+1/2 2 (3.8) 

where N is a normalizing constant. The diagonal operators A and C act on <]?£ as 

A<f> e = -l /2$t, h = m, Zi = 0, m = 0,1,2,... 

(3.9) 

C^ e = -h/2^i, h = n, n = 0,1,2,... 

This means that §i is a fundamental state of the representations j\ = m/2 of the subalgebra 
su(2) generated by {A^, A}, and j2 = n/2 of the corresponding su(2) determined by {C^, C}. 
Such a representation of su(3) will be denoted (m,n) with m,n G Z + . The points labeling 
the states of this representation obtained from lie on the plane D = m — n inside the 
^-parameter space. 

The energy for the states of the IUR's determined by the fundamental state (|3.9|l with 
the parameters (Zo, 0, Z2), according to (|3,6j) is given by 

E = (Z + h + 3/2)(Z + h + 5/2) = (m + n + 3/2) (m + n + 5/2) (3.10) 

Therefore, the IUR's fixed by (m, n) with the same value m + n will lead to states with the 
same energy. We call such IUR's an iso-energy series and they will be examined under the 
light of the algebra so(6) in the following section. The values for the energy (|3.1()j) coincide 
with the ones computed by the method of variable separation of section as can be seen 
from (|2.22j) once the replacement Zi = 1/2 is performed. We can also check that in this case 
the ground state (|2.24|) coincides with those fixing an IUR (|3,8|) . 

3.2 The so(6)— hierarchy 

Following the pattern and motivation of section 12.31 we will consider the relevant discrete 
symmetries in order to find a larger dynamical algebra. 

It is obvious that the Hamiltonian Hn 0: i 1: i 2 -\ is invariant under reflections in the parameter 
space {(Z , Zi, Z 2 )} 

Iq : (Zo, Zi , Z2) — > (— iQihih), h '■ {h-,h,h) — > {h,— h,h), h '■ {h,h,h) — > {h,h-,—h) 

Each of these symmetries can be directly implemented in the eigenfunction space, leading 
through conjugation to another set of intertwining operators that close a Lie algebra isomor- 
phic to u(3) and denoted by «u(3) 

j X = I i XI i , le«(3), iXe,u{3), i = 0,1,2 

The intertwining operators of iu(3) connect eigenstates of Hamiltonians whose parameters 
(Zq, Zi, Z2) belong to the planes {D = k{, being k{ certain real constants. We will choose the 
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following convention for the resulting generators 

{A ± ,B ± ,C ± } {A^B^C*} 

{A ± ,B ± ,C ± } {i ± ,S ± ,C ± } 

{A^S^C*}-^ {A^B^C^} 

where, for instance, the sets {A^A} and {A ± , A} close the two commuting Lie algebras 
su(2) of section |2J The explicit expression for the new operators (labelled with a tilde) can 
be easily obtained in the same way as it was done in (|2.18|) . The set of all the generators 
obtained in this process close the Lie algebra of rank 3, so(6). In the eigenfunction space it 
is enough to consider three independent diagonal operators {Lq, Li, L 2 } defined by 

The Hamiltonian can be expressed in terms of the so(6)-Casimir operator by means of the 
'symmetrization' of the u(3)-Hamiltonian ()3.6j) 

H so (6) = I ( H u(3) + Hj !j H u(3) Ij + T,j^k Ijh Su(3) Ijh + hhh H u( 2) I0W2) 
= {A+, A-} + B-} + {C+, C-} 

+ {A+,A-} + {B+, B-} + {C+, C-} + L 2 + ^i 2 + L 2 2 + f± 

Henceforth we remove the subindex 'so(6)' of the Hamiltonian. 

The intertwining generators of so(6) give rise to larger three-dimensional Hamiltonian 
hierarchies 

\H(la+m+p,h+m-n-p,l2+n)}i m,n,p € Z 

each one including a class of the previous ones coming from u(3). The eigenstates of these 
Hamiltonian hierarchies can be classified in terms of so(6)-representations. Let us fix our 
attention in those determined by the so(6) IUR's. These IUR's are build from the fundamental 
states annihilated by the simple roots A~ ,C~ , A~ 

A~*f = C-i,f = A~i,P = 

The equations for the operators A~ and C~ have been used in (|3.7|) . while the one for A~ 
were already applied in ()2.19|) . Therefore, the wavefunctions of the highest weight vectors 
take the form 

** 0) (</>i,02) = N cos 1/2 0i sin 1/2 0i cos02 sin i2+1/2 2 
characterised by the eigenvalues of the diagonal operators, 

L = L\ = 0, L 2 ^>i = n^i, n £ Z + 

This fundamental state is invariant under the inversions Iq and I±, and the representation, so 
obtained, is also invariant under I 2 . Thus, in this way we arrive at two classes of symmetric 
IUR's of so(6) 
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(a) (7 = 0, li = 0, l 2 = 0); {H {m+pm _ n _ pn) }, n,m,p eZ (even IUR's) 

(b) (l = 0, h = 0, l 2 = 1); {iJ( m+P)OT _ n _ Pi i +n )}, n,m,p G Z (odd IUR's) 

Each of these IUR's is described in the parameter space by an octahedral lattice of points 
such that it will include an iso-energy su{2>) (or jsu(3)) series of representations, quoted in 
the above subsection, which correspond to parallel exterior faces of the octahedron and some 
of its sections. Such sections are determined by the values of the diagonal operator D (or 
Di) whose values fix the corresponding u(3)-representations. 




Figure 1: Plot of the points representing the states of two odd IUR's with q = 1 (left) and 
q = 3 (right). The 6 (q = l)-eigenstates share the energy E = |-|. The 50 (q = 3)-eigenstates 
share the energy E = i ■ | (the points corresponding to q = 3 include those of q = 1, of the 
inner octahedron, which are doubly degenerated. 

For instance, the so(6)-representation labelled by n = 1, corresponding to the odd hier- 
archy, includes the first su(3)-series, (1,0) and (0,1) described by the opposite faces of an 
elemental octahedron. The so(6)-representation of the even hierarchy fixed by n = 2 includes 
the su(3)-series made of (2,0), (1,1), and (0,2). Those associated to (2,0) and (0,2) corre- 
spond to opposite triangular faces, while (1, 1) is described by the parallel hexagonal section 
through the origin. These features can be better appreciated in Figures ^ an d |U 

In general, the so(6) IUR's fixed by the parameter q will include the iso-energy series 
of the sw(3)-representations labelled by (m, n) with m + n = q. This is the degeneration 
explained by the larger algebra so(6). A similar discussion can be done with respect to the 
representations of the su(2) su(2) subalgebra. They can be identified with square sections 
of the octahedron. 

4 Eigenstates and factorizations 

Let Hi and be two Hamiltonians related by means of a differential operator X in the 
following form 

X H t = He X X ] E v = H e X ] (4.1) 
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Figure 2: The figure on the left is for the q = 1 IUR of so(6) where the triangular oposite 
faces correspond to two IUR's of su(3). The figure on the right correspond to the the points 
of a q = 3 IUR of so(6). The three sections describe three IUR's of su(3). 

where the dagger denotes adjoint differential operators. Then, it is said that X is an inter- 
twining operator connecting Hi with H^. 

In a formal way, the eigenfunctions of Hi are transformed by X into eigenfunctions of H^ , 
but one must be careful about the behaviour of some properties, such as square-integrability, 
singularities, or boundary conditions, which might be altered by X. The intertwining problem 
just as introduced in (|4.1|) . which applies to the u(3)-system of section |3 takes into account 
shape invariance, in the sense that the partner Hamiltonian differs from the initial Hi 
simply by changing the values of the parameters: In general, shape invariance leads 

to an algebraic structure of the intertwining operators as it happens in our present case. 

Now, we will discuss in this section the form of the su(3) intertwining operators of section|3] 
and its relation to certain eigenstates (similar considerations also apply to so(6)). First of 
all, note that we can write such operators (see expressions (j3.1|) . ()3.3j) and (|3.4j0 as 

Af = a ± + a l , Bf = b ± + (3 i , Cf = c ± + lt (4.2) 

where a^, b , (r 1 stand for vector fields (expressed, for instance, in the variables (f>\, fo) defined 
on the sphere and ai,Pi,~/i design functions also defined on the sphere. Notice that 

a+ = -(a")* = J 2 , b + = -(ft")* = h c+ = -(<f)t = J 

where Jo,Ji,J2 close the rotation algebra so(3). Moreover, taking the hermitian conjugate 
we have made use of the invariant measure on the sphere. If we write the Hamiltonians in 
the hierarchy displaying the kinetic (or free) part and the potential as 

Hi = # (kin) + Vi 
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we see that the vector fields originate the kinetic term, i.e. 

F (kin) = a + a - + b + b- + c+c 



and the components ai,(3i,^e (defined on the sphere) give rise to the potential Vf(0i,02), 
labelled by the parameters £ = (lo,li,h)- Substituting (j4,2j) in the Hamiltonian IjH.fiJ) and 
taking into account (j3.5|) . we get the expression 



where \i is a number depending on Iq,1\,12- Equation (|4,3[) can be considered as a non- 
linear partial differential equation linking the unknowns {ot£, fajji} with the potential, in a 
quite similar way to the Riccati equation for the superpotential uj in the one-dimensional 
Schrodinger equation. For this reason, we sometimes will refer to {otg,fa, 7^} as superpoten- 
tial functions. This is in agreement with a more general result where the first order 
intertwining is built by 'dressing' the symmetries of the Laplacian operator with certain 
functions. 

The basic property of the one-dimensional superpotential uj was that it could be consid- 
ered as the logarithmic derivative of a Hamiltonian eigenstate (see (|'2,1U[) ). Here, we have 
something similar with respect to the superpotentials {ae,fa,^i} but first we want to settle 
this problem in general terms. If we know an intertwining operator X satisfying (|4.1() it can 
help us in computing certain eigenfunctions of Hi. Notice that if we define the kernel, ICx, 
of X as the linear manifold of wave- functions annihilated by X, 



then, such a space is invariant under the Hamiltonian operator Hg. Thus, we can look for 
eigenfunctions inside fCx, in general a much simpler problem. But, in the case of X being 
a partial differential operator, its kernel includes certain arbitrary functions, so it is still an 
infinite dimensional space. This is in sharp contrast with ordinary first order differential 
operators where the kernel is one-dimensional. 

Another option we have at hand is the following. The intertwining relation (|4.1|) implies 
the commutation 



not necessarily that one annihilated by X, as was the case just considered above. In this 
case, however, a similar expression to (j'i.lOj) in terms of such eigenfunctions is no longer valid 
for uj. When we know several intertwining operators, as in the present case, we can apply 
them in different ways according to the above comments. 

i) Superpotentials associated to a global fundamental eigenstate of {A~ , B~ ,C~}. We 
consider the intersection of the kernels of all the intertwining operators. Assuming 
that this subspace is one-dimensional we have just one eigenstate (up to a factor) 
$0 annihilated by all the lowering operators {A~ , B~ , C~}. So that, we obtain the 
following expressions quite similar to (|2.10[) 



Viifafa) = (a?f + (a + a e ) + (fa) 2 + (b + fa) + ( le ) 2 + (c+ 7 ,) + A, 



(4.3) 



Xip = 0, 



X^XH,, = H t X ] X 

This means that we can look for eigenfunctions of Hi inside any eigenfunction space of X^X, 




(4.4) 



This mechanism corresponds to the IUR's characterized in section [3J 
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ii) Superpotentials associated to a partial fundamental eigenstate. If the above subspace 
is the trivial null space, we can still restrict ourselves to the kernel subspace of anyone 
of the intertwining operators, for example A~ . Thus, let <3? be an eigenfunction of Hg 
with 3> € /C^-, i.e. A~Q = 0. This allows us to set 

at = -(a"$)/$ 

From this equation we can also separate variables in <1>. So that, the eigenfunction 
equation H<& = E§ leads to a second order ordinary differential equation whose solution 
can be easily obtained. 

However, we must outline that in this case the remaining superpotential functions fy, j# 
have not a simultaneous expression (|4.4|) in terms of the same 3>, they need different 
eigenfunctions. Under this point of view, section [21 constitutes an illustration of how 
this option leads to eigenfunctions separated in the variables 4>i, fa- 
in) Other excited eigenstates. The second option is to solve, for instance, the eigenvalue 
problem A + A~& = a<&, requiring at the same time $ to be also a Hamiltonian eigen- 
function. In terms of the ambient coordinates so,si,S2 this equation is (see also 

The same procedure can be applied with other more general sets of operators commuting 
with the Hamiltonian. For instance, we can diagonalize H inside the subspace 

(e 2 A + A~ + ei B + B~ + e C + C~) $ = a $ 

where the e^'s are constant coefficients. This leads to eigenfunctions separated in elliptic 
coordinates, that we do not consider here 

5 Conclusions 

We have shown how to deal with the u{3) (and the general u(n) case (HI follows the same 
pattern |19j ) analog of a class of factorizable one-dimensional potentials with underlying dy- 
namical algebra u(2). The higher rank systems in consideration are well known inside the 
class of sup erintegr able Hamiltonians and, of course, our objective was not to compute orig- 
inal eigenfunctions. Our interest was to apply a different point of view to understand some 
properties in a new context. For instance, the classification of the irreducible representations 
of su(3) in series corresponding to so(6)-octahedrons, and the relations involved in this frame- 
work is a non trivial result that could be best appreciated inside the intertwining technique. 
The relation of the unitary representations with an special form of the superpotential func- 
tions, or the separable eigensolutions determined in terms of intertwining operators clarifies 
some of the known procedures. 

We have seen how the elements of one-dimensional factorizations must be adapted to the 
new context. For example, the relation of superpotentials and a whole class of eigenfunctions 
(not just one), the expression of the Hamiltonian operator is not just a simple factorization, 
the lattice of states must be drawn in a three-dimensional space, etc. 
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There are several problems that can be adressed using the present procedure. The systems 
underlying noncompact algebras u(p,q), inhomogeneous Lie algebras iu(p,q) and contracted 
algebras [2U] are among the first applications that we expect to report in a near future. But, 
in general, any other integrable Hamiltonian system will allow for this treatment, with or 
without variable separation. This application would be of most interest. 
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